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Theory  of  Field  Adsorption  of  Hydrogen 


Xiaoitiing  Ye  and  H.J.  Kreuzer 
Department  of  Physics,  Dalhousie  University 
Halifax,  N.S.  B3H  3J5,  Canada 


and 


D.R.  Salahub 

D'epanement  de  Chimie,  Universit'e  de  Montreal 
C.P.  6128,  Succ.  A,  Montreal,  Que.  H3C  3J7,  Canada 


Abstract 

Field  effects  on  the  adsorption  characteristics  of  H ,  and  on  metals  are 
calculated  using  density  functional  theory  for  a  jellium  metal  and  for  small 
metal  clusters.  New  dissociation  channels  are  explored  and  a  mechanism  for 
the  formation  of  a  field-stabilized  linear  species  is  proposed. 
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1.  Introduction 

Hydrogen  adsorbs  on  metal  and  semiconductor  surfaces  in  either  atomic  or 
molecular  form  in  on  top,  bridge  and  multicoordinated  hollow  sites.  Adsorp¬ 
tion  energies  for  atomic  hydrogen  range  from  2  to  4  eV;  molecular  hydrogen, 
on  the  other  hand  is  weakly  physisorbed  with  less  than  0.5  eV.O) 

Electric  fields  of  the  order  of  volts  per  angstrom,  i.e.  of  the  same  order 
of  magnitude  as  intra-atomic  fields,  dramatically  effect  the  electronic 
properties  of  the  surface  layers  of  metals  and  semiconductors For  in¬ 
stance,  on  clean  metals  electric  fields  reduce  the  activation  barrier  of 
evaporation,  eventually  to  zero  at  the  evaporation  field  strength.  Electric 
fields  also  change  the  binding  characteristics  of  adsorbates.  For  instance, 
the  binding  energies  of  rare  gases,  weakly  physisorbed  in  the  field-free 
case,  increase  by  at  least  an  order  of  magnitude  due  to  field-induced  chemi¬ 
sorption.  Likewise,  for  adsorbed  molecules  dissociation  may  become  enhanced 
or  inhibited  as  the  field  is  increased  changing  reaction  pathways.  Hydrogen 
is  no  exception  to  this  scenario.  Two  effects  in  particular  have  so  far 
defied  a  detailed  theoretical  explanation;  (/)  the  reduction  of  the  evapora¬ 
tion  field  strength  in  the  presence  of  an  electric  field,  and  (//)  the  appear¬ 
ance  of  a  species  in  the  field  ion  mass  spectrum^**- . 

In  this  paper,  we  will  report  theoretical  studies  of  the  adsorption  charac¬ 
teristics  of  atomic  and  molecular  hydrogen  as  a  function  of  electric  field 
strength.  Two  quite  general  field  effects  simplify  our  task:  (/)  Electric 
fields  are  always  enhanced  at  kink  sites,  terraces  and  in  front  of  single 
atoms  on  closed  packea  planes.  We  can  therefore  restrict  our  calculations 
to  on  top  sites.  (//)  Molecules  are  usually  alligned  along  the  field  direc¬ 
tion;  we  will  therefore  only  consider  such  geometries. 
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We  will  see  that  the  binding  energy  of  atomic  hydrogen,  H,  is  reduced  as 
the  field  is  increased,  whereas  is  bound  more  strongly.  As  a  result,  dis¬ 
sociative  adsorption  and  associative  desorption  are  modified  considerably. 
Within  our  theoretical  framework  we  will  also  be  able  to  propose  a  reaction 
pathway  for  the  formation  of  a  linear  species  at  metal  surfaces  in  high 
electric  fields.  A  detailed  analysis  of  experimental  data  in  the  light  of 
the  present  theoretical  results  will  be  given  elsewhere. 

2.  Theoretical  Models 

Electric  fields  at  metal  surfaces  are  greatly  enhanced  around  kink  sites  and 
around  single  metal  atoms  on  flat  crystal  planes.  These  local  electric 
fields  of  the  order  of  V /\  arise,  via  Poisson's  equation,  from  local  enhance¬ 
ments  of  the  surface  electron  density.  The  field  distribution  around  iso¬ 
lated  atoms  on  an  otherwise  flat  jellium  svurface  have  recently  been  calcu¬ 
lated  self  consistently  within  the  framework  of  the  local  density  approxima¬ 
tion  of  density  functional  theory A  field  enhancement  by  a  factor  1.5  to 
2  is  found  some  1-2%  above  an  adsorbed  metal  atom.  Calculated  evaporation 
field  strengths,  e.g.  2.QV/A  for  Ti  and  2.5V /A  for  Nb,  agree  very  well  with 
experimental  data. 

To  understand  field  adsorption  of  molecules,  one  would  like  to  calculate 
the  electronic  structure  of  an  adsorbed  molecule  in  the  local  field  at  kink 
sites  or  around  single  atoms  on  densely  packed  planes  of  transition  metals. 
Unfortunately,  such  a  complete  calculation  is  not  available  at  the  present 
time.  We  will  therefore  resort  to  two  approaches,  both  of  which  are  based  on 
density  functiop=>l  theory  to  determine  the  electronic  structure  and  the 
local  electric  field  selfccnsistently,  namely  the  jellium  model  to  describe 
field  adsorption  of  atomic  hydrogen  and  a  cluster  model  for  field  adsorption 
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of  the  mocelular  species. 

For  the  field  adsorption  of  atomic  hydrogen  we  will  use  the  jellium  model 
to  represent  a  densely  packed  metal  surface  on  which  we  will  study  the  chem¬ 
isorption  of  atomic  hydrogen  as  a  function  of  applied  field  using  local  den¬ 
sity  functional  theory  .(*2. i3)  one  assumes  that  the  ionic  lattice  can  be  smo¬ 
othed  into  a  uniform  positive  charge  density  /i+  that  drops  to  zero  abruptly 
half  a  lattice  constant  above  the  topmost  layer  of  ion  cores.  In  front  of 
this  jellivun  metal  we  add  a  hydrogen  atom^i**)  around  which  we  will  study  the 
local  electron  and  field  distribution  within  the  framework  of  density  func¬ 
tional  theory.  The  latter  determines  the  exact  ground-state  electron  den¬ 
sity  as  the  self-consistent  solution  of  a  Schrbdinger-like  equation  and  of 
Poisson's  equation.  The  local  density  approximation  is  invoked.'^)  To  in¬ 
clude  an  external  electric  field,  one  adds  an  excess  surface  charge  whose 
spatial  distribution  is  then  determined  self  consistently . 

To  account  for  the  local  lattice  effects  for  adsorption  around  isolated 
atoms  on  densely  packed  surfaces  or  at  kink  sites,  we  will,  in  a  second  set 
of  calculations,  use  a  finite  cluster  model  based  on  spin  density  functional 
theory.  Because  the  number  of  metal  atoms  is  rather  small  in  such  calcula¬ 
tions  for  practical  reasons,  long  range  metallic  effects  such  as  image 
forces  are  poorly  described  by  such  a  model.  The  inclusion  of  electric 
fields  is  also  no  longer  staightforward.  If  we  do  this  by  adding  charges  to 
the  cluster,  they  will  distribute  themselves  over  its  surface.  By  choosing 
an  appropriate  geometry,  e.g.  a  pyramidal  shape,  we  can  ensure  that  most  of 
the  extra  charge  resides  on  the  top  atom  of  the  cluster  and  a  strong  elec¬ 
tric  field  is  generated  in  its  vicinity.  This  field  will  fall  off  rather  rap¬ 
idly  like  r“2  a  few  angstroms  away,  rather  than  approach  a  constant  as  it 
does  in  front  of  a  single  metal  atom  on  top  of  an  extended  metal  surface. 


.  5  - 

Still,  one  can  generate  fields  of  about  IV /\  one  or  two  angstroms  away  from 
the  top  atom  in  the  cluster,  i.e.  at  the  position  where  atomic  hydrogen  will 
adsorb.  In  earlier  cluster  calculations<'5- based  on  the  semi-empirical 
ASED-MO  method(‘*)  the  electric  field  distribution  as  calculated  for  a  flat 
jellium  surfaces  was  imposed  on  the  cluster.  This  method  successfully 
explained  field- induced  chemisorption  of  rare  gases^'^-'"^  and  field  evapora¬ 
tion  of  metals^’®) .  For  the  cluster  calculations  reported  here  we  employ  a 
different,  and  simpler  scheme  to  include  field  effects.  Because  we  are  in¬ 
terested  in  field  adsorption  around  single  metal  atoms  on  top  of  closed 
packed  surfaces  or  around  kink  sites,  we  take  a  cluster  consisting  of  a 
single  metal  atom  and  the  adsorbing  atom  or  molecule.  This  cluster  is  then 
put  into  a  constant  external  field  which  is  supposed  to  arise  from  the  flat 
metal  surface  beneath  it.  The  induced  field  will  then  lead  to  a  net  field 
expulsion  over  the  volume  of  the  metal  atom  and  to  a  field  enhancement  in 
front  of  it  where  field  adsorption  of  the  additional  atoms  or  molecules 
takes  place. 

To  set  up  the  cluster  model,  we  again  start  from  the  one-electron  equa¬ 
tions  of  density  functional  theory  and  add  the  potential  energy  of  a  con¬ 
stant  electric  field,  Fq,  i.e.  (z)  =  -cFqZ  to  the  effective  potential.  For 
the  exchange  and  correlation  part  of  the  effective  potential  we  employ  both 
local  and  nonlocal  approximations.  To  allow  for  more  flexibility  we  divide 
the  total  electron  density  n{r)  into  spin-up  and  spin-down  components 

3 .  Results 

3.1.  H,  H2,  and  their  Ions 

To  assess  the  reliability  of  the  local  and  nonlocal  spin  density  approxima¬ 
tions  (LSDA  and  NLSDA),  we  have  compiled,  in  Table  1,  bond  lengths,  total 


-  6  - 


electronic  energies,  binding  energies,  and  ionization  energies  for  the  vari¬ 
ous  hydrogen  species  and  their  ions.  In  the  basis  set  for  hydrogen  we  in¬ 
cluded  Ls,  25,  and  2/?-type  gaussian  orbitals.  For  atomic  hydrogen  LSDA  is  not 
too  reliable  yielding  an  ionization  energy  of  13.02  eV  instead  of  the  exact 
value  of  13.59  eW',  a  similar  discrepancy  exists  for  the  ionization  energy  of 
//i-  the  molecular  species  bond  lengths  calculated  in  LSDA  are  typically 
too  long  by  as  much  as  0.1  Ionization  energies  for  linear  and  triangular 
Hj  are  in  good  agreement  with  Cl  calculations The  latter,  being  specif¬ 
ically  designed  for  ,  employ  much  larger  basis  sets  and  ought  to  be  more 
reliable.  We  also  note  that  Cl  calculations  based  on  the  Gauss' 90  package'--' 
give  rather  poor  results  as  does  an  approach  based  on  a  basis  of  floating 
spherical  Gaussian  orbitals .(^3)  Proceeding  to  the  NLSDA  we  note  a  signifi¬ 
cant  overall  improvement. 

As  a  first  step  towards  the  understanding  of  field  effects  on  hydrogen  we 
calculate  the  electronic  properties  of  in  a  constant  electric  field  using 
LSDA.  In  Fig.l  we  show  the  energy  difference,  AE-  E(H2,F)^2E{H,F=0)  as  a 
function  of  the  H-H  distance  for  various  field  strengths.  We  note  a  steady 
decrease  in  the  minimum  of  the  total  energy,  with  the  bond  length  increasing 
slightly.  The  Schottky  barrier  of  dissociation  disappears  at  5.14  k/H.  The 
earlier  FSGO  calculation  put  the  dissociation  field  strength  between  5.6  and 
6.0  To  assess  the  reliability  of  our  method  we  have  also  calculated 
the  field  dependence  of  the  activation  energy  of  dissociation  for  • 
According  to  Hiskes  <^4)  it  dissapears  for  fields  larger  than  2V /%.  whereas  we 
find  1.9V /%. 

3.2  Field  Adsorption  of  H 

For  the  field  adsorption  of  atomic  hydrogen  we  present  results  modeling  the 
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metal  by  (i)  a  jellium  and  (ii)  a  small  cluster  of  metal  atoms.  In  Fig.  2  we 
show  the  potential  energy  curves  of  H  on  jellivun  with  r^=2.Q  for  various 

field  strengths.  As  the  field  is  increased,  the  position  of  Lhe  binding  min- 

'  -^0 

imum  remains  essentially  unchanged  at  0.25A  above  the  jellium  edge,  i.e. 

0 

about  1  -  1.5  A  above  the  topmost  ion  cores.  The  activation  energy  for  field 
desorption  is,  however,  decreasing  continuously,  suggesting  that  atomic 
hydrogen  can  no  longer  be  field  adsorbed  at  field  strengths  beyond  5V/A. 
Slight  changes  occur  when  one  takes  r^=3.0  confirming  earlier  results  in  the 
field  free  case  which  found  that  the  total  energy  has  a  flat  plateau  for 
between  2  and  4  .(25.26) 

To  get  an  idea  for  the  importance  of  the  localized  p-  and  J-electrons  of 
transition  metals  and  of  the  lattice  structure  we  now  report  the  results  of 
LSDA  cluster  calculations  for  Rh.  In  a  first  calculation  we  have  mimicked 
the  (111)  surface  of  Rh  by  a  tetrahedral  cluster  of  4  Rh  atoms  with  a  H  atom 
on  its  apex.  It  is  known  that  atomic  hydrogen  adsorbs  preferentially  at 
multi-coordinated  sites  in  the  field-free  case.  Nevertheless  we  have  used 
the  on  top  geometry  because  in  the  presence  of  a  field,  adsorption  is 
preferrable  at  on  top  sites,  i.e.  at  kinks  and  seeps.  In  the  absence  of  a 
field  we  find  the  binding  energy  for  H  on  Rh^  to  be  2.36  eV  (compared  to  1.6eV 
on  jellium  with  r^=3  bohr  and  2.0  eV  for  r^=2  bohr)  with  H  located  1.5  ^  from 
the  top  Rh  atom  (compared  with  1.6^  on  jellium  if  we  assume  that  the  first 
lattice  plane  is  half  a  lattice  constant,  i.e.  l.'SrSA  below  the  jellium  edge). 
In  a  third  calculation,  also  based  on  LSDA,  we  look  at  the  Rh-H  dimer.  We 
find  a  binding  energy  of  3.3eV  at  a  distance  of  1.5  i.e.  somewhat  stronger 
binding  at  shorter  distances  for  the  dimer  than  for  the  cluster  because  the 
single  Rh  in  the  dimer  has  more  electrons  available  for  binding  to  hydrogen. 
The  further  reduction  in  binding  as  one  goes  to  jellium  is  largely  due  to  the 
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lack  of  localized  p-  and  £/-electrons . 

Within  the  dimer  model  we  next  include  field  effects  by  placing  the  dimer 
in  a  uniform  external  field.  This  overestimates  the  field  effect  on  Rh, 
because  of  the  absence  of  the  field  e:{pulsion  by  conduction  electrons  in  an 
extended  metal  cluster.  Nevertheless,  the  calculation  has  its  merits  as  an 
estimate.  In  Fig. 3  we  plot  the  potential  energy  curves  for  the  Rh-H  dimer  as 
a  function  of  the  separation  for  various  field  strengths.  It  suggests  that 
for  Fq  of  the  order  of  3. 0-3. 5 the  dimer  no  longer  exists.  This  is  indeed 
the  range  above  which  field  adsorption  of  hydrogen  on  metals  can  no  longer 
be  observed.  The  fact  that  a  dimer  model  yields  good  answers  confirms  a 
long  held  belief^-^-*^  that  the  local  fields  and  thus  the  local  electronic 
structure  around  adsorption  sites  (kink  sites,  steps  etc.)  are  more  impor¬ 
tant  than  long  range  metallic  effects  for  the  understanding  of  field  adsorp¬ 
tion.  It  is  at  first  sight  surprising  that  the  jellium  model  predicts  too 

large  maximum  field  strengths  for  hydrogen  field  adsorption,  i.e.  of  the 
0  0  .  . 

order  of  5-6V//1  instead  of  2-2. 5V/ A,  considering  the  fact  that  similar  cal¬ 
culations  for  field  evaporation  of  metal  atoms  yielded  good  agreement  with 
experiment .^1'^  However,  we  should  bear  in  mind  that  hydrogen  adsorption 
takes  place  at  very  short  distances  from  the  topmost  ion  core  of  the  metal, 
i.e.  in  the  region  of  the  classical  image  plane  where  the  local  field  varies 
rapidly.  The  local  structure  must  therefore  be  known  rather  acciirately,  a 
task  that  the  jellium  model  of  a  metal  cannot  provide.  The  geometries  and 
energies  of  field-free  adsorption  are  summarized  in  fig. 4  with  field  effects 
displayed  in  fig. 5. 

3.3  Field  adsorption  of 

We  have  used  the  LSDA  method  to  calculate  the  binding  energy  of  a  linear  Rh- 
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H-H  species  and  of  H2  linearly  adsorbed  on  top  of  a  tetrcihedron  of  Rh.  Such 
a  geometry  is  of  course  very  unlikely  in  the  absence  of  a  field.  However,  as 
soon  as  a  field  is  applied,  it  will  align  the  molecular  dipole  moment,  i.e. 
the  molecular  axis  along  it.  We  will  therefore  only  study  field  effects  for 
this  geometry. 

In  the  absence  of  a  field  we  find  a  binding  energy  in  LSDA  of  0.66  and 
0.54eV  on  a  single  Rh  atom  and  on  a  Rh^  tetrahedron  respectively,  i.e.  too 
much  binding  for  a  wecikly  physisorbed  species.  Using  NLSDA  the  binding  of  //, 
on  a  single  Rh  atom  decreases  to  a  more  realistic  value  of  0.065eV.  Unfortu¬ 
nately,  we  could  not  do  the  calculation  for  on  a  Rh^  cluster  due  to  rstic- 
tions  in  computer  size.  In  the  LSDA  calculations  the  first  hydrogen  is  a 
distance  of  1.88  (1.90)  %  from  the  top  Rh  atom  and  the  second  H  a  further  0.79 
%  away,  roughly  the  bonding  distance  in  the  free  molecule,  see  fig. 4. 
These  numbers  indicate  physisorption  for  molecular  hydrogen.  The  weakening 
of  the  Rh-‘H  bond  as  one  goes  form  a  single  Rh  to  a  Rh^  cluster  is  again  the 
result  of  a  small  charge  transfer  into  the  internal  Rh  bonds.  The  reduction 

of  the  binding  energy  as  one  employs  the  NLSDA  is  also  reflected  by  the  fact 

0 

that  the  binding  distance  increases  to  2.11  A. 

We  next  look  at  the  effect  of  a  uniform  electric  field  on  the  binding  char¬ 
acteristics  of  a  linear  Rh-H-H  species.  As  the  field  is  increased,  the  dis¬ 
tance  between  Rh  and  the  inner  H  shrinks  and  the  H-H  separation  increases 
enhancing  the  binding  between  the  Rh  and  the  H2  but  diminishing  the  H-H 
binding  energy.  For  fields  of  the  order  of  2V /%  the  activation  barrier  for 
H-,  dissociation  (in  a  direction  normal  to  the  surface)  dissappears.  Some 
relevant  numbers  can  be  found  in  fig. 5.  We  should  stress  that  this  is  a  new, 
field- induced  dissociation  channel  for  an  upright  //,  species  leading  to  only 
one  H  atom  adsorbed  rather  than  two  as  a  result  of  the  stretching  of  an  H^ 


-  10  - 


molecule  physisorbed  parallel  to  the  surface. 

3.4  Field  Adsorption  of 

The  bond  lengths  and  binding  energies  for  field-free  adsorption  of  a  linear 
molecule  on  either  a  Rh^  cluster  or  on  a  single  Rh  atom  are  given  in  fig. 4. 
This  suggests  to  interpret  adsorption  as  chemisorption  of  a  H  atom  on  Rh 
with  a  H2  molecule  physisorbed  on  top  of  it.  Dissociation  energies  support 
this,  i.e.  0.1  eV  for  Rh-H  plus  H2  and  3.28  eW  for  Rh  plus  Clearly  this  is 
not  the  lowest  energy  state  for  the  adsorption  of  3  hydrogen  atoms  on  a 
metal  surface. 

What  we  want  to  show  next  is  tha"'"  the  linear  R^  configuration  in  an  on  top 
position  will  be  stabilized  in  the  presence  of  an  electric  field.  We  there¬ 
fore  put  the  Rh-H-H-R  cluster  into  a  uniform  electric  field  of  2V/^  and  cal¬ 
culate  the  binding  energy  as  a  function  of  the  Rh-H^  distance  keeping  the  H- 
H  distance  at  its  field-free  value.  The  minimum  remains  at  the  same  dis¬ 
tance  but  the  activation  barrier  against  thermal  field  desorption  is  reduced 
to  less  than  leV.  If  we  vary  next  the  distance  between  the  first  and  the 
second  H  atoms  we  find  that  a  substantial  contraction  occurs  with  a  signifi¬ 
cant  increase  in  the  corresponding  H-H  bond  strength  implying  a  field  sta¬ 
bilisation  of  the  linear  molecule.  The  relevant  bond  lengths  and  binding 
energies  are  summarized  in  Figs.  4  and  5. 

4.  Summary 

The  main  conclusions  of  this  study  are  as  follows: 

(/)  Adsorption  of  atomic  hydrogen  on  metal  sxirfaces  is  wecikened  in  the  pres¬ 
ence  of  electric  fields.  The  maximum  field  strength  of  the  order  of  3  V/^  is 
reproduced  by  a  simple  dimer  model  whereas  field  effects  on  adsorbed 
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hydrogen  are  underestimated  on  a  jellium  surface. 

(//)  H2  linearly  adsorbed  in  an  on  top  position  will  be  dissociated  in  fields 
of  the  order  of  2V/\. 

(iii)  linearly  adsorbed  t.IH  be  stabilized  in  fields  of  the  order  of  2^/^. 

In  fields  of  less  than  3'//^  we  therefore  have  present  on  the  surface  (/) 
chemisorbed,  i.e.  immobile  atomic  hydrogen,  and  (//)  weakly,  field-stabilized 
chemisorbed  most  likely  mobile  even  at  temperatures  of  the  order  of  50K. 
The  H2  molecule  will  then  migrate  up  the  field  gradient  along  the  metal  tip. 
If  it  encounters  a  chemisorbed  H  atom,  it  might  hop  on  top  forming  a  field- 
stabilized  linear  species  which  can  be  field  desorbed  as  . 

We  hope  to  do  further  calculations  for  hydrogen  field  adsorption  in  three¬ 
fold  hollow  sites.  In  addition,  we  will  check  the  influence  of  the  metal 
cluster  size  on  the  field  adsorption  characteristics. 
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Table  1  Binding  characteristics  of  H,,  linear  and  triangular  and  their  ions. 
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Figure  Captions 


Fig.l; 


Fig. 2: 


Fig. 3: 


Fig. 4: 


Total  energy  of  H2  in  a  uniform  electric  field,  Fq,  as  a  function 
of  the  H-H  distance.  From  top  to  bottom,  the  field  strengths 
are  0.0,  3.0,  4.5,  5.14  Vj^A. 

Field  adsorption  of  atomic  hydrogen  on  a  jelliiam  metal  with 

r^=2.0.  From  top  to  bottom,  the  field  strengths  (at  the  minimum) 

0 

are  3.0,  5.0,  6.0,  0.0  VjA. 

Total  energy  of  the  Rh-H  dimer  in  a  constant  electric  field,  F^. 
From  top  to  bottom,  the  field  strengths  are  0.0,  2.0,  3.0,  3.5 
Vl%. 

Field-free  adsorption  of  hydrogen  on  jellium,  a  single  Rh  atom 
and  a  tetrahedron  Rh^.  Intemuclear  distances  are  in  angstroms 
and  energies  in  eV .  The  cluster  calculations  are  based  on  LSDA. 


Fig. 5: 


Field  adsorption  of  hydrogen  on  a  Rh  atom,  using  NLSDA. 
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